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Coherent control of quantum systems has far-reaching implications in q uant um engineering. In 
this context, coherent population trapping (CPT) involving dark resonanceP^ has played a promi¬ 
nent role, leading to a wealth of major ^plications including laser cooling of atom^and moleculeJ^, 
optical magnetometrjEl, light storagd^^® and highly precise atomic clockj^. Extending CPT meth¬ 
ods to individual solid-state quantum systems has been only achieved in cryogenic environments 
for electron spin impuritie^^^^ and superconducting circuit^^^^^^. Here, we demonstrate efficient 
CPT of a single nuclear spin in a room temperature solid. To this end, we make use of a three-level 
system with a A-configuration in the microwave domain, which consists of nuclear spin states ad¬ 
dressed through their hyperfine coupling to the electron spin of a single nitrogen-vacancy defect in 
diamond. Dark state pumping requires a relaxation mechanism which, in atomic systems, is sim¬ 
ply provided by spontaneous emission. In this work, the relaxation process is externally controlled 
through incoherent optical pumping and separated in time from consecutive coherent microwave 
excitations of the nuclear spin A-system. Such a pumping scheme with controlled relaxation allows 
us (i) to monitor the sequential accumulation of population into the dark state and (ii) to reach a 
new regime of CPT dynamics for which periodic arrays of dark resonances can be observed, owing to 
multiple constructive interferences. This work offers new prospects for quantum state preparation, 
information storage in hybrid quantum systemJ^SllII and metrology. 


Nuclear spins in solids have attracted considerable in¬ 
terest over the last years owing to their high level of iso¬ 
lation from the environment, leading to very long coher¬ 
ence timePHU. 

The detection and control of individual 
nuclear spins is usually achieved by exploiting their hy¬ 
perfine coupling to an ancillary electronic spin. This can 
be perfo rmed by using a nitrogen-vacancy (NV) defect in 
diamoncPHHl whose spin triplet ground state {S = 1) can 
be optically initialized, coherently manipulated with mi¬ 
crowave helds and read-out by optical means through its 
spin-dependent photoluminescence (PL) intensity^. Ex¬ 
periments making use of the hyperfine coupling between 
the NV defect and nearby long-lived nuclear spins are 
now reaching a level of control allowing to realize elabo¬ 
rate quantum information protocol m°”^^l Here, we make 
use of electron spin transitions of a single NV defect to 
achieve coherent control of a single nuclear spin through 
CPT under ambient conditions. 

The spin system considered in this study is depicted in 
Fig.[T^. The electronic spin of a single NV defect is cou¬ 
pled by hyperfine interaction with a nearby nuclear 
spin (/ = 1/2). In the rris = 0 (|0e)) electron spin mani¬ 
fold, the nuclear spin eigenstates are denoted |t) and j),), 
corresponding to the nuclear spin projections along the 
NV defect quantization axis. In the TOs = —1 (| —le)) 
electron spin manifold, efficient nuclear spin mixing can 
be obtained by exploiting the anisotropic component of 
the hyperhne interaction in combination with a magnetic 
field applied along the NV defect axid^. The nuclear spin 




FIG. 1: Spin system, a, A single NV defect is coupled with 
a nuclear spin via hyperfine interaction. The energy level 
structure shows the nuclear spin eigenstates in the ms = 0 
(|0e)) and ms = —1 (| —E)) electron spin manifolds. The 
arrows indicate electron spin resonance (ESR) transitions, b, 
ESR spectrum of the spin system recorded at a magnetic field 
B — 850 G by sweeping the frequency of a microwave field 
while monitoring the spin-dependent PL intensity of the NV 
defect. Here the nuclear spin mixing parameter is 0 « 1.14 rad 
(see Supplementary Information). 
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FIG. 2: Coherent population trapping of a single nuclear spin, a, A-scheme configuration. Relaxation from the excited 
state is triggered by optical pumping with a green laser, b, CPT pulse sequence. Two selective MW rr-pulses of duration 
Tmw = 6 /rs are followed by a 300-ns-long green laser pulse, with a power corresponding to the saturation of the NV defect 
optical transition. These parameters ensure spectral selectivity of an isolated A-scheme and enable optimizing the electronic 
spin readout contraslP^. c, CPT dip observed in the ESR spectrum while setting di = 0. For each sweep, the CPT sequence is 
repeated N ~ 2000 times before incrementing <52. d, ESR spectra recorded for various two-photon resonance conditions. The 
red solid lines are data fitting with two inverted Gaussian functions. All experiments are performed under ambient conditions 
with a magnetic field B — 850 G applied along the NV axis. 


eigenstates are then given by 

1 + ) = cos(^0|t)+sin(^^)e*^|i) 

I-) = -sin +COS li) , 

where the mixing parameter 9 can be tuned with the 
magnetic field strength and the phase </> is fixed by the 
hyperfine tensor components [see Supplementary Infor¬ 
mation]. In this case, four transitions can be observed in 
the electron spin resonance (ESR) spectrum, with rela¬ 
tive intensities fixed by the nuclear spin mixing parame- 
teiUll [Fig.[^]. The spin system therefore exhibits two A- 
level configurations in the microwave domain which can 
be used to perform CPT of the nuclear spin. 

As illustrated in Fig. one of these A-schemes was 
isolated and its ESR transitions coherently excited with 
two microwave fields, MWI and MW2, detuned from the 
resonances by (5i and 62 - To briefly remind the general 
principle of CPT, we consider the two orthogonal linear 
superpositions of the lower levels 


where Hi and H 2 are the Rabi frequencies characteriz¬ 
ing the strength of the MW coupling on each branch of 
the A-scheme. At the two-photon resonance condition 
<5 r = (5i — (^2 = 0, the population of the dark state \D) is 
fully decoupled from the MW fields. On the other hand, 
the population of the bright state |H) is efficiently cou¬ 
pled to the excited state, and then partially transferred to 
the dark state through spontaneous decay. After few ex¬ 
citation cycles, the system is thus trapped into the dark 
state. Such a CPT process requires an efficient relax¬ 
ation from the excited state. For A-schemes in the opti¬ 
cal domain, this is provided by the irreversible process of 
spontaneous emissiorP. In this work, the intrinsic decay 
from the excited state | —) is governed by the longitudi¬ 
nal electron spin relaxation of the NV defect. However, 
this process occurs within a typical timescale Ti of few 
milliseconds at room temperature and is bi-directional, 
thus leading to a thermal state. To circumvent these 
detrimental effects, we use optical pumping to polarize 
the NV defect electronic spin in jOe). More precisely, the 
population transfer to the dark state is deterministically 
triggered within a sub-microsecond timescale by apply¬ 
ing a green laser pulse. This is the key ingredient of our 
work. 

The experimental sequence used to observe CPT of a 
single nuclear spin is depicted in Fig. Two MW pulses 
driving selectively the two branches of the A-scheme are 
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followed by a green laser pulse, which is used both for 
pumping the spin system into the dark state and for op¬ 
tical readout of the electronic spin state. ESR spectra 
were then recorded by continuously repeating this CPT 
sequence while monitoring the spin dependent PL inten¬ 
sity of the NV defect and sweeping the detuning 82 across 
the resonance, with (5i being fixed. A typical ESR spec¬ 
trum recorded for = 0 is shown in Fig. revealing 
a pronounced CPT dip when the two-photon resonance 
condition is fulfilled. This is the signature of an efficient 
pumping of the spin system into the dark state. Re¬ 
markably, the contrast of the CPT dip is close to unity 
(« 90%), which implies that the laser pulse used to trig¬ 
ger the relaxation process is not altering significantly the 
nuclear spin coherence. Repeating the experiment for 
different values of shows that the position of the CPT 
dip shifts in accordance with the two-photon resonance 
condition = 0 [Fig. [^], confirming the dark state 
pumping interpretation. 

Such a pulsed CPT scheme gives a unique opportunity 
to analyze the sequential accumulation of population into 
the dark stat^^ while increasing the number N of rep¬ 
etitions of the CPT sequence. To this end, the two mi¬ 
crowave fields were applied resonantly (Ji = ^2 = 0) and 
the probability V\d){N) to find the spin system in the 
dark state was inferred by measuring the population in 
state | —) at each repetition of the CPT sequence through 
a calibrated spin-dependent PL intensity measurement 
(see Supplementary Informations). The experimental re¬ 
sults are shown in Fig. |^. Starting from a thermal state 
^|D> (0) = 1/2, the dark state population increases expo¬ 
nentially with the number of CPT steps and saturates to 
V\D){'X)) ~ 0.88 ± 0.03, which corresponds to the dark 
state polarization efficiency. This dynamics is well re¬ 
produced by a simple rate equation model describing the 
competition between dark state pumping with a proba¬ 
bility Qfp per step and depolarization induced by nuclear 
spin decoherence with a probability a^p per step. The 
probability ap mainly depends on the MW pulse areeP^ 
and the unbalanced relaxation between the two branches 
of the A-scheme, while the probability adp is essentially 
linked to nuclear spin decoherence induced by optical 
pumping (see Supplementary Information). Here, data 
fitting leads to ap = 0.43 ± 0.07, adp = 0.12 ± 0.03 and 
the characteristic number of steps Ng of the exponential 
growth is A^s = 1.4 ± 0.3. The dark state composition 
at the steady state was further analyzed by measuring 
1(1 |H)P, i.e. the probability to find the spin system in 
state ID, as a function of the Rabi frequencies (Hi, 112) 
of the microwave fields driving the A-system. As shown 
in Fig. 1^, this probability scales as VL\/{n\ -|- Hf), as 
expected from the dark state expression in the {|t), ID} 
basis [see Eq. &■ This is a further evidence of CPT of 
the nuclear spin. 

As already pointed out, our pulsed-CPT experiment 
is performed in an unusual regime where the A-scheme 
relaxation is (i) deterministically triggered and (ii) well 
separated in time from the coherent excitation. These 




H1/H2 


FIG. 3: Step-by-step dark state pumping, a, Probability 
VlD) to find the nuclear spin in the dark state as a function 
of the number N of CPT steps. The experiment is performed 
with Hi « H 2 using resonant microwave 7r-pulses. The solid 
line is data fitting with a rate equation model, leading to a 
characteristic number of steps Ng of the exponential growth 
Ns = 1.4 ± 0.3. b, Analysis of the dark state composition 
at the steady state {N = 20) as the function of the Rabi 
frequencies ratio r = H 1 /H 2 . The solid line is data fitting 
with the function r^/(l -b r^) weighted by the imperfect CPT 
contrast (see Supplementary Information). 


properties enable to precisely tune the time Tgeq between 
subsequent coherent MW excitations [Fig. |^] and ana¬ 
lyze its effect on the properties of the dark resonance. 
For that purpose, ESR spectra similar to the one shown 
in Fig. were recorded while increasing the duration 
of the CPT sequence. An array of CPT dips appears 
in the spectrum, whose period A evolves as = Tgeq 
[Fig-i^ ,d]. Moreover, the width of the CPT resonances 
gets significantly narrower when Tgeq increases [Fig. He]. 
To understand these observations, we first analyze Ime 
CPT pulse sequence in the spectral domain. As shown 
in Fig. 1^, the periodically repeated MW pulses of dura¬ 
tion Tmw can be described by frequency combs, whose 
peaks are separated by T~^ within a carrier-envelope 
of spectral width In addition, the width of each 

peak of the comb is given by [AgTseq]”^, where Ng is the 
characteristic number of CPT steps needed to reach the 
steady state of the spin system. Using this simple fre¬ 
quency comb picture, it is straightforward to show that 
dark state pumping is possible whenever the two-photon 















4 


T 

seq 


MW1 


repeat 






FIG. 4: Observation of multiple dark-resonances, a, The total duration of the CPT sequence Taeq is modified by varying 
the waiting time between coherent MW excitations and relaxation induced by optical pumping, b, ESR spectra recorded in 
the same way as in Fig. n while increasing Tgeq (from top to bottom). The red solid lines are Gaussian fits used to extract the 
periodicity and the width of the GPT dips, c, Pulsed-MW excitations described as frequency combs driving the A-scheme. d, 
Frequency of the dark resonance array as a function of Taeq. The red solid line corresponds to = Taeq. e, Full width 
at half maximum (FWHM) of the central GPT dip as a function of Taeq. The red solid line is data fitting with the function 
[AsTaeq]”^, leading to W = 1.8 ± 0.1. 


detuning is such that (5 r, = nT^gq, where n is an inte¬ 
ger [Fig. 0 = 1 - The spectral interval A between consecu¬ 
tive dark resonances then follows = Tgeq, as exper¬ 
imentally observed [Fig. &]■ Furthermore, the width of 
the CPT resonances closely follows the one of the comb 
peaks given by [AgTseq]”^ Data fitting leads to 

Ng = 1.8 ± 0 . 1 , in fair agreement with the step-by-step 
measurements reported in Fig. 

These multiple dark resonances can also be understood 
by considering the evolution of the dark state in the time- 
domain. For a two-photon detuning Jr, it evolves as 

\m ) = (^2 It) - It)) , (3) 

corresponding to a precession at the frequency Jr be¬ 
tween the states \D) and |i?) defined by Eq. ([^. When 
the phase accumulated within Tseq is such that Jr = 
nTggq, the dark state is locked in phase with the peri¬ 
odic CPT pulse sequence, leading to efficient dark state 
pumping through multiple constructive interference^^. 

The dynamics of pulsed CPT with controlled relax¬ 
ation markedly contrasts with continuous (CW) dark 
state pumping schemes. Indeed, CPT spectra recorded 
in CW regime usually result from a competition between 


the excited state relaxation rate 7 and the coherent cou¬ 
pling. This is manifested by the fact that the CPT width 
is given by D 2/7 in the limit where <C ^2 and neglect- 
ing decoherence of the ground state superpositioiP^l. In 
this work, the relaxation is temporally decoupled from 
coherent coupling, and is not playing a role in the width 
of the resonance. In this novel regime, the width of the 
CPT dip is rather fixed by the sequence duration and 
Ns, which could be efficiently tuned by changing the 
area of the microwave pulses, as theoretically analyzed 
in Ref. n. 

The fundamental limit to the CPT width is determined 
by two intrinsic parameters of the spin system. First, 
dark state resonances can be observed as long as Tgeq 
is shorter than the longitudinal electron spin relaxation 
time Ti, which is on the order of few ms at room tem¬ 
perature, and can be extended to several tens of seconds 
at cryogenic temperatur^^. This limits the CPT width 
to Fi = [AgTi]”^. Moreover, Tgeq needs to be shorter 
than the dephasing time of the nuclear spin = F^^, 
which can be as long as few seconds under ambient con¬ 
dition^^. The fundamental limit to the CPT linewidth 
is therefore given by Fg = max{ri,r 2 „}, which could 
potentially reach the sub-Hz regime at low temperature. 
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In conclusion, we have demonstrated CPT of a sin¬ 
gle nuclear spin in a room temperature solid. The 
experiment is realized in a unique scenario for which 
dark state pumping is performed in a pulsed fashion 
while controlling the relaxation of a A-scheme in the mi¬ 
crowave domain. This novel CPT regime may lead to 
new approaches for coherent control of individual nuclear 
spins, e.g. through stimulated Raman adiabiatic passage 


(STIRAP), infor matio n storage and processing in hybrid 
quantum system^^MS] precision atomic clock^ 

We thank E. Arimondo, P. Pillet, D. Clement, P. 
Bertet, P. Maletinsky and J. R. Maze for fruitful dis¬ 
cussions. This research has been partially funded by the 
French National Research Agency (ANR) through the 
projects Advice, Qinvc and Smequi. 


SUPPLEMENTARY INFORMATION 


I. EXPERIMENTAL SETUP 

We study native NV defects hosted in a commercial [100]-oriented high purity diamond crystal grown by chemical 
vapor deposition (Elementd) with a natural abundance of ^^C isotopes (1.1%). Individual NV defects are optically 
isolated at room temperature using a home-built scanning confocal microscope under green laser excitation. Coherent 
manipulation of the NV defect electron spin is performed by applying a microwave field through a copper microwire 
directly spanned on the diamond surface. A permanent magnet mounted on a xyz-translation stage is used to apply 
a static magnetic field along the NV defect axis. Details about the experimental setup can be found in Ref. [5^ . 


II. SPIN HAMILTONIAN 


The NV defect in diamond has a spin triplet ground state S' = 1, whose spin projection along the NV quantization 
axis (z) is denoted as rus = 0, ±1. In this section, we do not introduce the hyperfine interaction between the NV defect 
electronic spin and its intrinsic ^^N nuclear spin (/jv = 1 ), since it is not playing any role in the reported experiments. 
We therefore consider a fixed ^“^N nuclear spin projection, e.g. mj^.^ = -|-1. For a magnetic field amplitude B applied 
along the NV defect axis, the ground-state spin Hamiltonian of the NV defect in frequency units reads 

no = DSl + , (4) 


where D « 2.87 GHz is the zero-field splitting and 7 e ~ 2.8 MHz/G is the electronic spin gyromagnetic ratio. When 
a neighbouring lattice site of the NV defect is occupied by a isotope, corresponding to a nuclear spin / = 1/2, 
the spin Hamiltonian becomes 


n = no + inBi, + s • A • i, 


(5) 


where 7 „ = 1.07 kHz/G is the gyromagnetic ratio of the ^^G nuclear spin and A its hyperfine tensor. In the secular 
approximation, it simplifies as 

n = Ho + InBI, + + e*^/_) , (6) 

where Aani = {Aj^ + tancj) = A^y/Azx and /± = 4 ± ily. 

In the following we consider that Azz > 0, and the eigenstates of the 4 operator are denoted as |t) and H). As 
indicated in the main text, the eigenstates ji/i) of the spin system in the TOs = 0 (|0e)) and nia = —1 (| — Ig)) electronic 
spin manifolds are given bjl^ 


IV'l) = | 0 e,t) 

|■^/' 2 ) = |0e,4.) 

(^0 A^I-leA) 

11/4) = I - le,-) =-sin - le,t)+COS 


IV'a) = 


- le, -k) = cos ( - ) I - le, t) + sin 


(7) 

( 8 ) 

(9) 


l-leA) 


( 10 ) 
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FIG. 5: (a)-Energy level structure of the spin system showing the nuclear spin eigenstates in the = 0 (|0e)) and rria = ±1 
(I ± le)) electron spin manifolds. We consider Azz ~ 7nB and Azz > 0. The colored arrows indicate electron spin transitions 
labelled from ® to (6). (b)-ESR spectrum recorded for the coupled spin system with a magnetic field B = 850 G applied along 
the NV defect axis. The low frequency part of the spectrum is reproduced from Fig. 1 of the main text. 


with 


The corresponding eigenenergies read 


tan 9 = 




Azz -InB 


InB 

InB 

l/o = - 

2 

^3 = D--feB- + {Azz - InBy 

Vi = D-'yeB+ + (Azz - InBY . 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


When Azz InB, i.e. for 9 ~ 7r/2, four transitions from rus = 0 to rrig = — 1 can be induced by a microwave field, 
owing to nuclear spin mixing in the iris = —1 electronic spin manifold [see Fig. [^. Within this four level system, a 
A-scheme can be isolated, e.g. {|0e,t)j |0e,4.), | — le, —)}, in order to perform CPT of the nuclear spin . 

For all the experiments reported in the main text, we use a nearby nuclear spin whose hyperfine interaction is 
characterized by a longitudinal component Azz ~ 1 MHz and an anisotropic component Aani = {A?^x + A^^y^A ~ 
0.3 MHz. It corresponds to the lattice site denoted as I in Ref. [M]. A magnetic field B = 850 G is applied along 
the NV axis, which ensures efficient nuclear spin mixing in the mg = — I electronic spin manifold, while remaining far 
away from the ground-state level anti-crosssing {B — 1024 G). At this field, the nuclear spin mixing is characterized 
by 0 « 1.14 rad. The coherence time of the NV defect used in this work is « 8 /iS at H = 850 G, which is long 
enough to allow selective excitation of the A-system. 
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In the rris = +1 (|+le)) electron spin manifold, which is not considered in the main text, the eigenstates of the 
spin system are given by 


IV's) — I + le) +0 — cos ( —^ I + le,|) + sin + le, t) 


IV'e) = I + le,-') =-sin ( — ) e *'^| + le,4.) + cos 


+ le, t) I 


with 


tan0' = 


A.a 


•^zz + InB 


Here {Azz + 7 ™^) > Aani, so that 0' « 0 and 


with the corresponding eigenenergies 


IV's) = I + le4) 

IV'e) = I + le,t) ) 


1^5 = D +-feB - -{Azz + InB) 


1^6 — D + JeB + — {Azz + 7 ™-®) ■ 


(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 


The nuclear spin quantization axis is therefore parallel to the NV defect axis and only two nuclear-spin con¬ 
serving transitions from = 0 to mg = -1-1 can be observed in the ESR spectrum [Fig. §>]■ These high-frequency 
ESR transitions were used to calibrate the experiments of sequential dark state pumping as explained in Section[lII C| 


III. SEQUENTIAL DARK STATE PUMPING 
A. Rate equation model 

In this section, we analyze the experimental results reported in Figure 3 of the main text by introducing a simple rate 
equation model, which describes how the populations evolve within the A system while repeating the CPT sequence. 
We consider N successive dark state pumping steps with two microwave fields applied at resonance, i5i = (52 = 0 [Fig. 
[^a)]. The spin system is described in the basis {|T1), \B), |—)}, where \D) and \B) are the dark and bright states of 
the A system, which are defined by Eqs. (1) and (2) in the main text. Each pumping step can be decomposed into 
three different parts, as shown in Fig. |^b). 

1. The transition \B) —>■ |—) is first coherently driven though pulsed MW excitation with an effective Rabi frequency 

Q = a pulse area A = Hr, where r is the MW pulse duration. 

2. The spin system is then kept in the dark during a waiting time shorter than the longitudinal relaxation time Ti 
of the electronic spin. Here, the populations of the A system do not evolve as long as the nuclear spin dephasing 
in the dark can be neglected. 

3. A laser pulse is then applied at t = Tiasen which triggers the relaxation from the excited state |—) to the ground 
states, \D) and \B), with rates T/j and Tb respectively. At the same time, the laser pulse induces nuclear spin 
dephasing with a rate Tdp. The laser pulse duration is denoted as AT = Tggq — Tjaser- 

We first isolate the pumping step number N. The populations of the ground states at the beginning of step N are 
equal to the ones obtained at the end of step {N — 1), which are labelled V\d) — 1) and V\b) — 1) for state \D) 
and |H), respectively. Considering an ideal polarization of the NV defect electron spin under optical illumination, 
these populations are linked by the relation V\d){^ — 1) + 'B\b){N — 1) = 1. After coherent MW excitation, the 
population in state |—) is given by 


T|_) {N) = sin^ V^B) {N-l) = sin^ [l - T|b> {N - 1)] 


(23) 







FIG. 6: (a)-Pulse sequence used for dark state pumping. (b)-Evolution of the populations in the A-system (represented as 
orange spheres) during the CPT sequence. Q, = is the effective Rabi frequency on the transition \B) —>■ |—) , Td and 

r B are the optically-induced decay rates and P^p is the dephasing rate of the nuclear spin under optical illumination. 


During optical illumination, dark state pumping competes with laser-induced dephasing of the nuclear spin with a 
rate We denote by ap the probability to decay from |—) to the dark state at each CPT step. This probability is 
linked to the decay rates of the A system through 

= rap = r|(-|D)|2 (24) 

Fb = ^(l-ap) = ^|(-|S)|^ (25) 

where F is the optically-induced polarization rate of the NV center electronic spin in the TOs = 0 stat^^. 

The time evolution of the populations during the laser pulse, i.e. for t S [Tjaser, Tseq], is then described by the rate 
equations 


dVl-) 

dt 

d'P\D) 

dt 

dViB) 

dt 


r^i-> 

-^{'PiD)-riB)) + rapVi.^ 
^{VlD)-V^B)) + r{l-ap)Vi.^ . 


Considering a closed A-system also implies that 


(26) 

(27) 

(28) 


P\D) (t) + P\B) (t) + P\-) (t) = 1. 

By solving equation (26) for t G [Tiaser, Tseq] with the initial condition at t = Tiaser given by Eq.(23), one can find 

lP|_)(f) = lP|-)(riaser)e-^(‘-^'“-^ = sin^ [l - IPp) (IV - 1)] ). 


Using Eqs. (29) and ( |30[ ), the evolution of the dark state population becomes 
dP\D) 


dt 


leading to 


P\D){t) — 


+ TdpP\D) (t) = P\-) (Tla«er)e-^(‘-^'—) 
P\-) (Tiaser) 


Foip 


F 


dp 


Fdp 


F-F 


dp 


r 0!.r) 


dp 


g — rdp(i —Tiaser) g— T(i— Tiaser) 


Tid) (Tiaser) ^ 


„ —rdp(t —Tiaeer) 


(29) 
ind 

(30) 


(31) 


1 

2 


(32) 
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where 'P\d) {T\ aser )=V\d){N-1). 

At the end of the laser pulse of duration AT = Tgeq — Tiaser, we consider that the electronic spin is fully polarized 
into ms = 0, so that « 0. The link between the dark state population at the end of the CPT sequence 

'^\D) (Tseq) = 'P\D) [N) and the one at the beginning of the step 'P\d) {N — 1) then reads as 


ViD)iN)=ViD){N-l)e 


— Fdp AT 


1 - 


sin^ (I 




+ 1 ri _ 

2 L ^ T-Tdp I 


(33) 




dp 


By recurrence and with an initial thermal state of the nuclear spin [T|£)>(0) =0.5] 
N steps is finally expressed as 


, the dark state population after 


'Pid){N) 


b 


1 — a 


H” d 


N 


7^|d)(0) 



with 


a 


b 


g PdpAT ^ 


T — Tdp 




1 [1 _ e 


T — r^p 




(34) 


(35) 

(36) 


The characteristic number of steps needed to reach the stationary value of the dark state population is given by 


Ns = -l/ln(a) = -1/ln 


1 - 


sin' 


hi) 


r-r, 


dp 


Fa - hp 
lap 2 


(37) 


and the dark state population at the steady state reads 


'P\D)ioo) = 


1 [1 _ e-r.pAT] + 


F—Fdc 


Pft — 

1 Up 2 


X _ g-FdpAT 


si„^(4) 

^ r-r^p 


Fn- — 

^U'p 2 


(38) 


In order to simplify these equations, one can consider that dark state pumping is faster than laser-induced dephasing, 
i.e. QfpT ^ Tdp. In our experiment, this assumption is supported by the observation of a large contrast of the CPT 
dip in the ESR spectrum, which is the signature of an efficient pumping of the spin system into the dark state [Fig. 
2 of the main text]. Hence, the dark state population after N pumping steps simplifies as 


-Pid^N) 


a'p + adp{0.5 - a'p) ^^(I - a^p) N\ 

a'p + adp{l - a'p) ' -|-adp(l - op V Ns) 


(39) 


where a), = OpSin^ (Di ^dp = [1 ~ exp(—T^pAT)] is the depolarization probability per step induced by nuclear spin 
decoherence and Ns = —1/ln [l — a/ — adp(l ~ ci/)] is the characteristic number of pumping steps needed to reach 
the steady state of the A system. At the steady state, the dark state population is then given by 


P\D){oo) 


a'p + adp{0-5 - a'p) 
a'p + adpil - a'p) 


(40) 


corresponding to a regime for which the gain by a pumping step into the dark state equals the losses induced by 
depolarization. If a/ 3> adp, then 'P|d)(oo) = 1 and the system is fully polarized into the dark state. On the other 
hand, if adp ^ a/, then V\d){oo) = 0.5 and the system remains in a thermal state. Hence, effective dark state 
pumping can be achieved as long as a'p > adp- By fitting the step-by-step CPT experiments presented in Fig. 3 of the 
main text with equation (39), we obtain a'p = 0.43 (0.07) and adp = 0.12 (0.03), leading to V\D){ao) = 0-88 (0.03). 
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B. Polarization efficiency 


We now describe more precisely the polarization process. For clarity purpose, we assume that adp = 0 and that the 
MW pulses bring the whole population of the bright state \B) to the excited state |—), i.e. A = tt. We consider the 
asymmetric A system represented in Figure [^b). As a reminder, the spin states of the A-scheme are expressed by 


I-) 

1^) 

\B) 


-sm I - I e 
1 

y/ni + 

1 

y'ni + nj 


^'^|t)+cos(^0U) 


(41) 

(42) 

(43) 


where fli and 142 are the Rabi frequencies of the two MWs driving the A-scheme. Here we have introduced the phase 
difference ^ between the two MW fields, which was omitted in the main text. 

|Il)p corresponding to the probability to decay from 


The polarization efficiency is defined in Eq.(24) as Op = 


|—) to the dark state during laser illumination. IJsing Eqs. (41) and (42), one can obtain 


(Xn — 


1 -|- (141/122) 


sin X + 


(14i/142)2 
1 + (14i/142)2 


cos^ X + 


(14i/142) 

1 -|- (14i/142)' 


■ sin(d) cos{(j) — ■*/) 


(44) 


The parameters 14i, 142 and tp allow for a full control over the dark state polarization efficiency and the A-scheme 
relaxation [see Fig. [^. This might enable to modify in a controlled fashion the dynamics of step-by-step dark state 
pumping. 


(a) e = ■Kjl 


(b) 0 = 1.14 rad 




FIG. 7: Polarization coefficient calculated from Equation (44) for (a) Q = -k 12 and (b) Q = 1.14 rad corresponding to the nuclear 
spin mixing parameter for a magnetic field B = 850 G applied along the NV axis [see Section]^. In each case, the polarization 
coefficient is inferred for three different values of A4> = p — ip (A4> = 0 in blue, A4> = 7r/2 in black and A<1> = tt in red). 


C. Experimental procedure 

In this section, we give details of the experimental procedure which enables to extract the probability 'P\d){N) 
through a calibrated spin dependent PL intensity measurement. 

Reference sequence - The high-frequency ESR transitions, presented at the end of Section |TTj were used to 
calibrate the experiments of sequential dark state pumping and to measure the nuclear spin’s populations in [f) 
and [/). Assuming that the PL intensity is identical for TOs = -1-1 or TOs = —1, we use the sequence depicted on 
Fig. [^a) to make a correspondence between PL intensity and population. The NV defect electronic spin state 
is first polarized in nig = 0 by laser illumination. After a 1-^s long waiting time, a 300-ns laser pulse is applied 
and the recorded PL signal is used as a reference signal for the nig = 0 spin sublevel. A selective MW tt pulse 
is then applied on the transition [Oe,).) —>■ | -I- le,i) [MWhf in Fig. [^], which brings half of the population into 
rUg = +1. A 300-ns long laser pulse is then applied and the recorded PL signal is used as a reference signal for the 
nig = ±1 spin sublevels. These measurements were used to calibrate the experiments of sequential dark state pumping. 
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FIG. 8: (a)-Sequence used to make the correspondence between PL intensity and populations of the spin system. The microwave 
frequency MWhf is set on resonance with the transition |0e,4-) I + le,4') [labeled ® on figurej^c)]. (b)-Experimental pulse 
sequence used for step-by-step dark state pumping. 


Step-by-step sequence - As described in the previous paragraph, the experiment starts by applying a reference 
sequence which enables making the correspondence between PL intensity and populations of the spin system. Dark 
state pumping steps are then applied [Fig. [^b)]. Each pumping step is composed by two 6-/rs-long MW pulses 
driving the A-scheme at the two-photon resonance condition. Here the MW power is adjusted in order to obtain 
TT-pulses, as verified by recording electron spin Rabi oscillations. After a 100-ns-long waiting time, a 300-ns laser pulse 
is used to pump the spin system into the dark state. The laser pulse is followed by a 1-^s waiting time which ensures 
efficient relaxation of the population from the metastable state of the NV defect. The pumping steps are repeated 
N times and the population in state |—) of the A-scheme, i.e. in rUg = —1, is measured at each step by recording 
the spin-dependent PL during the laser pulses. After N repetitions, a 100 fxs laser pulse is used to reinitialize the 
system into a thermal state. The whole sequence is repeated continuously in order to improve the signal to noise ratio. 


Extracting the dark state population - The PL intensity measurement is done during each laser pulse of the 
CPT sequence. The recorded PL is compared to the reference signal to convert it in a measurement of the excited 
state |—) population, which follows Eq. (23). Considering an initial thermal state of the nuclear spin [V\d){0) — 0.5], 
we have 




(45) 


Using Eqs. (23) and (45), the 


probability to find the spin system in the dark state then follows 


Vid){N-1) = 1 


27 ^ 1 -) (!)■ 


(46) 


The sequential accumulation of population in the dark state can therefore be inferred by simply measuring the 
population in state |—), i.e. in rus = —1, at each pumping step through a calibrated spin-dependent PL intensity 
measurement. 


D. Characterization of the dark state 

The experimental sequence used to characterize the dark state is sketched in Figure The dark state’s composition 
is tuned by changing the power, and thus the Rabi frequencies (Di and D 2 ), of the MW fields driving the A system. 
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Laser _ 

MWhf - 

MWl - 
MW2 ^ 


Reference 

PL 


20 


6 us 






Dark state pumping 


TT 


2/is 


Population 

measurement 


300ns 

300ns 

100)18 






relaxation 


FIG. 9: Pulse sequence used for characterizing the dark state composition at the steady state. 


The pumping sequence is repeated TV = 20 times in order to reach the steady state population of the dark state. 
The projection of the nuclear spin along the NV center axis is then measured by applying a selective 7r-pulse on 
the high-frequency ESR transition |0e,i) —>■ | -h le,i) [MWhf in Fig. [^]. The spin-dependent PL is measured by a 
300 ns laser pulse and gives the probability to find the nuclear spin in state jj,). A 100 /is-long laser pulse is finally 
applied in order to relax the system to a thermal state. The whole sequence is repeated continuously in order to 
improve the signal to noise ratio. 

This measurement gives the probability to find the nuclear spin in state jj,). By representing the nuclear spin state 
with a density matrix pnuc, the measured probability (j, \pnuc\ i) at the steady state is expressed as 

{i\Pnuc\i) = V\D){^m\D)\^ + ii-ViD)i^)m\B)\^ . (47) 


Hence, one can find 


la |D)I’ = . (48) 

2V\d){oo) - 1 

In Figure 3(b) of the main paper, we plot |(), \D)\'^ as a function of H 1 /H 2 . Ideally, when the system is in the dark 
state, the probability to find the nuclear spin in state j),) is given by 

To take into account the imperfections of the measurement, the experimental data were fitted with the function 

/(Hi/H2)=0.5 + a(^^^2^-0-5) , (50) 

where a is an artificial contrast of the curve. Data fitting leads to a = 0.78 ± 0.03. 
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